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Abstract 


This  technical  memorandum  presents  two  improved  boundary  element  methods  for  pre¬ 
dicting  hydrodynamic  coefficients  of  two  dimensional  sections.  The  first  method  is  a  source 
distribution  method,  while  the  second  is  a  direct  method  for  determining  velocity  potential.  A 
major  objective  of  these  improved  methods  is  to  eliminate  errors  due  to  irregular  frequencies, 
which  currently  exist  in  the  close-fit  method  used  by  the  strip  theory  program  SHIPMO.  Both 
new  methods  successfully  eliminate  problems  caused  by  horizontal  and  vertical  line  segments 
with  the  close-fit  method.  To  eliminate  irregular  frequencies  of  surface-piercing  sections,  the 
source  distribution  method  employs  additional  sources  on  a  deck  lid  at  the  waterline.  The 
introduction  of  a  deck  lid  eliminates  some  irregular  frequencies  but  introduces  new  irregular 
frequencies.  Alternatively,  a  modified  direct  method  effectively  eliminates  all  irregular  frequen¬ 
cies.  Future  strip  theory  implementations  should  incorporate  this  modified  direct  method.  This 
technical  memorandum  also  reviews  low  and  high  frequency  linoits  of  hydrodynamic  coefficients. 

Resume 

Cette  note  technique  presente  deux  methodes  ameliorees  qui  resolvent  des  problemes  aux 
frontieres  pour  predire  les  coefficients  hydrodynamiques  de  sections  bidimensionnelles.  La  pre¬ 
miere  methode  utilise  une  distribution  de  sources  tandis  que  la  seconde  est  une  methode  di- 
recte  qui  determine  le  potentiel  des  vitesses.  Un  objectif  tres  important  pour  ameliorer  ces 
methodes  consiste  a  eliminer  les  erreurs  dues  aux  frequences  irregulieres,  qui  sont  actuellement 
presentes  dans  la  methode  d’ajustement  que  le  programme  SHIPMO  utilise.  Les  deux  nouvelles 
methodes  eliminent  avec  succes  les  problemes  que  la  methode  d’ajustement  rencontre  avec  les 
segments  de  lignes  horizontaux  ou  verticaux.  Pour  eliminer  les  frequences  irregulieres  dues  aux 
sections  qui  percent  la  surface  libre  la  methode  par  distribution  de  sources  utilise  des  sources 
supplementaires  disposees  sur  un  pont  au  niveau  de  la  ligne  de  flottaison.  Cette  technique  elimine 
des  frequences  irregulieres  mais  en  fait  apparaitre  d’autres.  Une  autre  solution  consiste  a  utiliser 
une  methode  directe  modifiee  qui  elimine  reellement  les  frequences  irregulieres.  Les  prochaines 
implementations  de  la  methode  de  decompoition  devraient  comprendre  cette  methode  directe 
modifiee.  Cette  note  technique  montre  egalement  les  limites  basses  et  hautes  frequences  des 
coefficients  hydrodynamiques. 
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EXECUTIVE  SUMMARY 


Introduction 

The  strip  theory  code  SHIPMO  is  DND’s  primary  tool  for  predicting  motions  and  loads  of 
ships  in  waves.  The  slow  development  of  alternative  tools  such  as  time  domain  codes  suggests 
that  strip  theory  will  continue  to  be  the  primary  seakeeping  tool  for  the  next  several  years.  One  of 
the  greatest  problems  for  strip  theory  computations  is  the  presence  of  irregular  frequencies,  which 
can  cause  very  large  errors  in  hydrodynamic  coefficients  for  surface-piercing  sections.  These 
errors  adversely  affect  motion  predictions  and  can  severely  affect  load  predictions.  This  study 
examines  two  improved  boundary  element  methods,  the  source  distribution  method  and  the 
direct  method,  and  their  abilities  for  removing  irregular  frequencies.  Prediction  of  hydrodynamic 
coefficients  at  the  extremes  of  the  frequency  range  is  also  examined. 


Principal  Results 

The  improved  source  distribution  and  direct  methods  eliminate  many  problems  with  the 
current  close-fit  method.  The  source  distribution  method  with  a  deck  lid  removes  some  irregular 
frequencies  but  introduces  new  irregular  frequencies.  In  contrast,  a  modified  direct  method  effec¬ 
tively  removes  all  irregular  frequencies.  At  high  frequencies,  the  infinite  frequency  limit  provides 
an  excellent  approximation  for  velocity  potential  and  added  mass.  Two  dimensional  methods 
give  unrealistically  high  heave  added  mass  and  heave  damping  values  for  surface-piercing  ship 
sections  at  low  frequencies.  A  suitably  chosen  cut-off  frequency  appears  to  provide  an  acceptable 
approximation  of  heave  added  mass  at  lower  frequencies. 


Significance  of  Results 

The  modified  direct  method  provides  more  reliable  predictions  of  hydrodynamic  coefficients 
than  the  close-fit  method  currently  used  in  SHIPMO.  Its  implementation  into  SHIPMO  will  give 
improved  ship  motion  and  load  predictions.  The  infinite  frequency  limit  provides  reliable  values 
at  high  ship  frequencies  and  can  also  be  applied  to  impulse  phenomena  such  as  ship  slamming. 
Better  estimates  of  low  frequency  hydrodynamic  coefficients  will  lead  to  better  estimates  of  ship 
motions  and  loads  at  low  encounter  frequencies. 


Ill 


Future  Plans 


The  modified  direct  method  will  replace  the  close-fit  method  in  the  next  version  of  SHIPMO. 
The  infinite  frequency  limit  for  hydrodynamic  coefficients  at  high  encounter  frequency  and  the 
improved  low  frequency  predictions  will  also  likely  be  implemented. 
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Notation 


[A] 
B 

[B] 

[C] 


{D} 

[E] 

Ei{z,) 

Ei{Z) 
Fpv{Y,Z) 
G{x,Xs) 
G{x,Xs) 
G{x^Xs) 
Go(x,  Xs) 
®s) 

9 

Kzz 

k 

L 

Ik 

Is 

n 

Uy^  Tlz 

n(^) 


'^max 


n 

Sh 

T 

V 

Vn 

w,- 


matrix  relating  velocity  potential  to  source  strength 
sectional  breadth  at  waterline 
matrix  relating  normal  velocity  to  source  strength 
matrix  relating  velocity  potential  to  vector  {D} 
wave  velocity 

velocity  term  vector  for  direct  method 
Green  function  matrix  for  direct  method 
exponential  integral  function 
exponential  integral  function 

principal  value  function  component  of  Green  function 
Green  function 

frequency  independent  component  of  Green  function 
frequency  dependent  component  of  Green  function 
Green  function  at  zero  frequency  limit 
Green  function  at  infinite  frequency  limit 
gravitational  acceleration 

coefficient  for  evaluating  heave  added  mass  of  ellipsoid 

wavenumber  or  mode 

ship  (or  ellipsoid)  length 

length  of  line  segment  k 

source  length 

normal  pointing  into  body 
normal  pointing  into  body  at  x 
y  and  z  components  of  normal  n 
normal  component  for  mode  k 
radius 

maximum  distance  from  roll  origin  to  sectional  offset  points 

radius  from  field  point  to  image  source 

body  surface 

sectional  draft 

velocity 

normal  velocity 

normal  velocity  for  mode  k  (time  independent) 
complex  normal  velocity  for  mode  k  (time  independent) 
Gaussian  weighting  factor  for  point  j 
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Notation  (continued) 

X 

field  point  or  source  point 

Xi,X2 

end  points  of  field  line  segment 

X* 

lateral  image  point  of  x 

Xs 

source  point 

Xsl’)  Xs2 

end  points  of  source  line  segment 

Xs-j 

quadrature  point  j  on  source  line  segment 

Y,Z 

dimensionless  coordinates  of  field  point  relative  to  image  source 

y,z 

horizontal  and  vertical  coordinates 

Umax 

maximum  sectional  y  value 

complex  variable  -Z  +  iY 

^9 

height  of  roll  axis  above  waterline 

a 

weight  term  for  modified  direct  method 

o^i 

angle  of  field  point  i  relative  to  source  point 

Oisi 

angle  of  source  point  i  relative  to  field  point 

A 

ship  (or  ellipsoid)  displacement 

Sij 

Kroenecker  delta  function 

C 

water  surface  elevation 

V 

body  displacement 

Vk 

body  displacement  for  mode  k 

0 

field  line  inclination  angle 

9s 

source  line  inclination  angle 

^ij 

damping  coefficient 

y- 

Green  function  integration  variable 

IJ'ij 

added  mass  coefficient 

i 

field  point  or  source  point 

p 

water  density  ' 

a 

source  strength  j 

(t(^) 

source  strength  for  mode  k  j 

$ 

velocity  potential 

4> 

complex  velocity  potential  (time  independent) 

complex  velocity  potential  for  mode  k 

<!>' 

revised  velocity  potential  for  low  and  high  frequency  limits 

u> 

oscillation  frequency 

(jJq 

low  frequency  limit  for  valid  heave  coefiicients 

(4) 

^0  ^ 

lowest  frequency  for  roll  coefficient  computations 

viii 
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1  Introduction 


Strip  theory  ship  motion  programs  such  as  SHIPMO  [1,  2,  3, 4,  5,  6]  are  the  most  commonly 
used  method  for  predicting  ship  motions  and  sea  loads  in  a  seaway.  A  key  aspect  of  a  strip  theory 
program  is  the  prediction  of  hydrodynamic  coefficients  of  two  dimensional  sections.  SHIPMO 
can  predict  sectional  hydrodynamic  coefficients  using  either  the  conformal  mapping  method  or 
the  Frank  close-fit  method  [7,  8]. 

The  conformal  mapping  method  has  various  practical  limitations.  The  geometry  of  an 
input  section  is  customarily  described  by  beam,  draft,  and  area  coefficient.  Consequently,  the 
method  can  only  model  sections  which  pierce  the  free  surface  and  fall  within  an  acceptable  range 
of  these  input  parameters.  The  limitation  of  using  only  three  geometric  parameters  is  that  the 
conformal  mapping  method  gives  a  limited  approximation  to  the  actual  sectional  geometry.  This 
approximation  can  significantly  affect  roll  added  mass  and  damping  predictions,  which  are  very 
sensitive  to  sectional  geometry.  Another  limitation  of  the  conformal  mapping  method  is  that  it 
does  not  predict  sectional  velocity  potentials;  thus,  long  wave  approximations  must  be  used  for 
predicting  sectional  excitation  forces  during  strip  theory  calculations.  However,  an  important 
advantage  of  the  conformal  mapping  method  is  its  robustness. 

Geometrical  input  for  the  Frank  close-fit  method  is  given  in  the  form  of  sectional  offsets, 
permitting  high  accuracy  in  input  geometry.  The  Frank  close-fit  method  used  by  SHIPMO 
suffers  from  problems  as  described  in  References  6  and  9.  Input  offsets  for  a  section  must  have 
successively  greater  values  for  the  vertical  coordinate  z\  thus,  the  code  is  not  able  to  handle  a 
section  with  a  flat  bottom.  The  required  spacing  between  vertical  coordinates  appears  to  depend 
partly  on  which  computer  the  code  is  run,  indicating  a  lack  of  robustness  in  the  algorithm.  A 
major  problem  of  the  Frank  close-fit  method  is  the  presence  of  irregular  frequencies  at  which 
results  are  very  erratic.  Irregular  frequencies  are  common  among  source  distribution  methods 
for  predicting  hydrodynamic  coefficients  of  a  surface-piercing  body.  The  irregular  frequencies 
are  caused  by  numerical  solutions  which  satisfy  prescribed  boundary  conditions  but  initiate 
sloshing  within  the  body.  This  phenomenon  occurs  at  the  upper  frequency  range  of  interest  for 
ship  motions. 

This  technical  memorandum  examines  two  improved  boundary  element  methods  for  pre¬ 
dicting  sectional  hydrodynamic  coefficients.  The  first  method  is  a  modified  source  distribution 
method  in  which  a  line  source  is  placed  across  the  waterline  of  a  section  to  inhibit  sloshing.  Such 
methods  have  been  applied  previously  by  Bedel  [10],  Ohmatsu  [11],  and  Ando  [12].  The  second 
method  is  a  direct  method  as  described  by  Sclavounos  and  Lee  [13].  The  method  of  Reference  13 
is  implemented  in  the  DREA  submarine  motion  code  SUBMO  [14,  15]  without  correction  for 
irregular  frequencies;  the  version  described  here  includes  the  correction. 


2  General  Description  of  Problem 

Strip  theory  requires  a  solution  to  the  flow  in  the  vicinity  of  an  oscillating  two  dimensional 
section  such  as  that  shown  in  Figure  1.  The  problem  of  the  flow  around  a  general  oscillating 
body  is  described  in  many  references,  including  Kim  [16],  Sarpkaya  and  Isaacson  [17],  and  Mei 
[18].  The  problem  is  formulated  in  terms  of  an  oscillatory  velocity  potential  $.  For  the  two 


dimensional  case,  motions  are  limited  to  sway,  heave,  and  roll.  The  body  motions  are  taken  as 
being  about  the  center  of  gravity,  located  a  distance  Zg  above  the  waterline. 


Figure  1:  Two  Dimensional  Ship  Section 


The  velocity  potential  $  must  satisfy  a  continuity  equation  and  several  boundary  equations. 
Laplace’s  continuity  equation  is: 


V2$  =  0 


(2.1) 


Bernoulli’s  equation  for  conservation  of  energy  must  be  satisfied  in  the  fluid.  Retaining  only 
linear  terms,  conservation  of  energy  imposes  the  following  boundary  condition  on  the  free  surface: 

■^+5C  =  0  at  2  =  0  (2.2) 

where  C  is  the  elevation  of  waves  radiated  from  the  body.  The  relationship  between  wave 
elevation  and  the  velocity  potential  on  the  surface  is: 


dz  dt 


at  2  =  0 


(2.3) 


From  Equations  (2.2)  and  (2.3),  the  velocity  potential  on  the  free  surface  must  satisfy: 


df^  ^  ^  dz 


at  2  =  0 


(2.4) 


Imposing  Equations  (2.2),  (2.3),  and  (2.4)  at  2  =  0,  the  undisturbed  surface,  is  a  consequence  of 
linearization.  The  following  radiation  condition  is  based  on  the  requirement  that  waves  generated 
by  the  oscillating  body  must  radiate  outward  from  the  body: 


1 

dr  ^  c  dt 


=  0 


at  r  =  c» 


(2.5) 


where  r  is  distance  from  the  body.  The  vertical  water  velocity  must  be  zero  at  the  water  bottom. 
For  deep  water,  this  condition  is: 


d^  - 

=  0  at  2  =  —00 

dz 


(2.6) 
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on  Sb 


The  final  boundary  condition  is  that  the  flow  into  the  body  must  be  zero: 

II  =  K  on  5-6  (2.7) 

where  n  is  the  normal  pointing  into  the  body  and  Vn  is  the  normal  velocity  of  the  body. 

For  a  body  experiencing  oscillatory  motions  at  frequency  w,  solution  of  the  problem  is 
greatly  simplified  using  a  time  independent  complex  velocity  potential  4>  which  has  the  following 
relationship  with  the  time  dependent  potential  #: 

=  Real{^(2/,  z)  expiut}  (2.8) 

Thus,  the  complex  velocity  potential  (f)  is  merely  a  function  of  space. 

In  a  manner  similar  to  that  used  for  the  velocity  potential,  the  normal  velocity  of  the  body 
can  be  expressed  as  a  complex  quantity  which  is  independent  of  time: 

V{t)  =  Real{tr  expioJt}  (2.9) 

where  V  and  v  are  the  actual  and  time  independent  oscillatory  velocities.  For  a  given  motion 
mode  k,  the  relationship  between  the  complex  velocity  and  the  complex  displacement 
is: 


(2.10) 


To  obtain  hydrodynamic  coefficients  for  a  two  dimensional  section,  the  complex  velocity 
potential  ^  around  the  body  must  be  solved  for  sway,  heave,  and  roll  (denoted  by  modes  k  =  2, 
3,  and  4  respectively).  From  Equations  (2.7),  (2.8),  and  (2.9),  the  complex  potentials  must 
satisfy  the  following  body  boundary  conditions: 


^ on  Sb  for  A:  =  2, 3, 4 

On 

(2.11) 

The  complex  normal  velocity  component  Vn^  for  mode  k  is  related  to  the  body  normal  component 

7j(*)  by: 

(2.12) 

For  motions  about  the  center  of  gravity,  the  normals  for  each  mode  are: 

=  ny 

(2.13) 

=  Uz 

(2.14) 

=  y  Uz  {z-  Zg)  Uy 

(2.15) 

Once  the  velocity  potentials  on  the  body  are  known,  the  added  mass 

and  damping  are 

easily  evaluated  from: 

fiij  =  -  1  Imag|<^(-'^j  Ui  dS 

(2.16) 

Xij  =  p  J  Real{<^(^)}  ni  dS 

(2.17) 

where  and  A,j  are  added  mass  and  damping  coefficients  for  mode  i  due  to  motion  of  mode  j. 


3 


3  Source  Distribution  Method 


Source  distribution  methods  for  predicting  two  dimensional  hydrodynamic  coefficients  are 
described  by  many  authors,  including  Kim  [16],  Sarpkaya  and  Isaacson  [17],  and  Mei  [18].  The 
Frank  close-fit  method  [7,  8]  is  an  example  of  a  source  distribution  method. 

The  flow  in  a  fluid  domain  can  be  described  using  a  distribution  of  complex  sources.  For 
flow  induced  by  a  two  dimensional  oscillating  body,  the  following  equation  can  be  used: 


where  x  is  location  in  the  fluid  domain,  Xs  is  source  location,  G{x,  Xs)  is  the  Green  function 
describing  the  flow  at  x  caused  by  a  source  of  unit  strength  at  x*,  and  a{xs)  is  the  strength  of 
the  source  at  Xg- 

The  frequency  domain  Green  function  for  a  two  dimensional  source  based  on  Newman  [19] 
is: 


G{x,Xs)  =  ln(r)  -  In(r-i) 

/•OO 

—  2  ®  - -  exp  [— /iZ]  cos(/iy)  dfi  -b  i  27r  exp{—Z)  cosY  (3-2) 

Jo  -  1 

The  sign  of  the  imaginary  term  in  the  above  Green  function  is  opposite  to  that  given  by  Newman 
but  is  consistent  with  Kim  [16]  and  Salvesen,  Tuck,  and  Faltinsen  [20].  The  terms  r  and  ri 
respectively  denote  the  distance  from  the  field  point  x  to  the  source  at  Xs  and  to  the  image  of 
the  source  above  the  free  surface: 

r  =  y/{y-  Vsf  +  (z-  Zsf  (3.3) 

ri  =  y/iy-Vsf  +  {z  +  Zsf  (3.4) 

The  dimensionless  coordinates  Y  and  Z  are  based  on  wavenumber,  are  both  positive,  and  are 
defined  as  follows: 

A;  =  —  (3.5) 

9 

Y  =  k\y-y,\  (3.6) 

Z  =  -k{z  +  zs)  (3.7) 

A  key  feature  of  the  Green  function  of  Equation  (3.2)  is  that  it  satisfies  the  required  conditions 
of  Equations  (2.1),  (2.4),  (2.5),  and  (2.6).  The  only  remaining  condition  is  the  velocity  boundary 
condition  on  the  body  (Equation  (2.11)),  which  can  be  satisfied  using  the  correct  source  strengths 
cr{x)  on  the  body  surface.  The  equation  for  determining  the  source  strengths  that  satisfy  the 
body  velocity  boundary  condition  is: 


4 


3.1  Discretized  Solution  of  Source  Distribution  Method 


Equation  (3.8)  is  normally  solved  by  discretizing  a  ship  section  into  several  line  segments, 
as  shown  in  Figure  2.  Ten  line  segments  are  usually  sufficient  to  describe  one  half  of  a  typical 
ship  section.  The  source  strength  a  is  assumed  to  be  constant  over  each  line  segment  and  is 
such  that  the  velocity  boundary  condition  of  Equation  (2.11)  is  satisfied  at  the  midpoint  of  each 
line  segment.  The  resulting  discretized  form  of  Equation  (3.8)  is: 


-  <T(fi)  + 


2N 

j=l 


dG{xi,Xs) 

dn^. 


(t{Xs)  dS  =  Vn{Xi) 


(3.9) 


where  2N  is  the  number  of  segments  into  which  the  ship  section  has  been  discretized,  Sj  is  the 
surface  of  segment  j,  and  x,-  is  the  midpoint  of  segment  i.  Note  that  the  number  of  segments  in 
Equation  (3.9)  is  denoted  2N  rather  than  N,  which  will  simplify  the  development  for  symmetrical 
sections  presented  later.  Once  the  source  strengths  cr(xs)  are  known,  the  discretized  velocity 
potentials  can  be  solved  using: 

^  /  G{xi,Xs)  (7{xs)  dS  (3.10) 


Figure  2:  Discretized  Ship  Section  for  Source  Distribution  Method 


The  matrix  forms  of  the  equations  relating  the  velocity  potentials  and  normal  velocities  to 
the  source  strengths  are: 


{#)}  =  [A]  {cr«} 
=  [B]  {ctW} 

From  Equations  (3.9)  and  (3.10),  the  matrix  coefficients  are: 

^ij  ~  o  f  a::^)  dS 

27r  Js^ 


B. 


=  - + -  / 
•'  2  ^  2ir  Is 


dG{^x^ ,  X5) 

dns: 


dS 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


where  6ij  is  the  Kroenecker  delta  function. 
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3.2  Symmetry  Relations  for  Source  Distribution  Method 


When  running  a  linear  ship  motion  program  such  as  SHIPMO,  lateral  symmetry  usually 
exists  for  ship  sections,  giving  the  following  ; 


4‘)(r)  = 

- 

for  A:  =  2, 4 

(3.15) 

vP(x-)  = 

for  A;  =  3 

(3.16) 

d>(^)(r)  = 

for  A;  =  2, 4 

(3.17) 

<f>^^\x) 

for  A:  =  3 

(3.18) 

<r('')(r)  = 

-  C7(*)(X) 

for  A;  =  2, 4 

(3.19) 

a('=)(r)  = 

(7(^)(X) 

for  A:  =  3 

(3.20) 

where  x*  is  the  lateral  image  of  x  defined  by: 

f  =  -V  (3.21) 

z*  =  2  (3.22) 

Using  the  symmetry  relations,  the  source  strengths  only  need  to  be  solved  on  one  side  of  the 
body.  If  the  source  strengths  and  velocity  potentials  are  solved  for  y  >  0,  then  the  coefficients 
for  Equation  (3.11)  and  Equation  (3.12)  are: 


(k) 


A'- 


Ak) 


51" 


(fc) 


Z3 


for  A;  =  2, 4 


for  A;  =  3 


dG{xi,Xs) 

dnsi 


dG{xi^  Xs) 


dug. 


dS  - 

dS  + 


dG{xi,x^) 

dnsi 

dG{xi,x:) 

dug. 


dS 

dS 


for  A;  =  2, 4 


for  A:  =  3 


(3.23) 

(3.24) 

(3.25) 

(3.26) 


where  SJ  denotes  the  lateral  image  of  Sj. 

Once  the  velocity  potentials  are  known  for  a  given  half-section,  the  added  mass  and  damping 
are  easily  computed  using  the  following  discretized  forms  of  Equations  (2.16)  and  (2.17): 


N 


flij 

^5  £  4Imag|4^^}  4'^ 

A;=l 

(3.27) 

Xij  = 

2p  AfcReall^-'^}  4^ 

k=l 

(3.28) 

where  Ik  is  the  length  of  segment  k. 
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3.3  Correction  of  Irregular  Frequencies 

A  major  problem  with  source  distribution  methods  is  the  occurrence  of  irregular  frequen¬ 
cies  at  which  Equation  (3.1)  gives  incorrect  values  for  source  strengths  on  an  oscillating  body. 
Irregular  frequencies  are  caused  by  sloshing  modes  which  can  occur  inside  a  surface-piercing 
body  because  the  formulation  (Equations  2.1  to  2.7)  allows  for  an  internal  free  surface.  Frank 
[7]  gives  the  following  equation  for  the  irregular  frequencies  of  a  heaving  rectangular  body: 


u>j  =  Y  coth(j  -K  T IB)  for  j  =  1,2, . . .,  oo  (3.29) 

For  sectional  computations  using  SHIPMO,  the  lowest  irregular  frequencies  (j  =  1,2)  typically 
occur  at  the  upper  portion  of  the  frequency  range  of  interest;  thus,  it  is  very  desirable  from  a 
practical  viewpoint  to  eliminate  irregular  frequency  effects. 

A  common  method  for  eliminating  sloshing  and  its  associated  irregular  frequencies  is  to 
place  a  deck  lid  across  the  waterline  of  a  surface- piercing  body,  as  described  by  Bedel  [10], 
Ohmatsu  [11],  and  Ando  [12].  Figure  3  shows  a  ship  section  with  a  deck  lid.  The  deck  lid 
introduces  additional  line  sources.  An  arbitrary  boundary  condition  can  be  imposed  on  the  deck 
lid.  Generally,  this  boundary  condition  must  be  an  odd  function  of  y  for  sway  and  roll  and 
an  even  function  of  y  for  heave.  In  practice,  the  boundary  condition  on  the  deck  lid  is  usually 
specified  to  be  that  the  normal  velocity  is  zero.  A  variable  number  of  sources  can  be  used  to 
model  the  deck  lid.  Computations  by  Bedel  [10]  and  Ando  [12]  indicate  that  a  small  number  of 
line  segments  (one  or  two)  on  the  -fj/  portion  of  the  deck  lid  gives  the  best  results. 


Figure  3:  Discretized  Ship  Section  with  Deck  Lid 


4  Direct  Method 

The  direct  method  is  an  alternative  to  the  source  distribution  method,  and  has  been  de¬ 
scribed  by  Sclavounos  and  Lee  [13]  and  Sarpkaya  and  Isaacson  [17].  Like  the  source  distribution 
method,  the  direct  method  has  problems  with  irregular  frequencies;  however,  Sclavounos  and 
Lee  indicate  that  the  errors  caused  by  irregular  frequencies  are  less  severe. 
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The  direct  method  permits  the  direct  solution  of  velocity  potentials  as  follows: 

..,dGi(,x)  _  If  dm 


i  di  =  —  f 

dn.  ^  2t:  Js 


iSt  dn^ 


(4.1) 


where  ^  is  a  point  on  the  body.  To  permit  solution  of  Equation  (4.1),  the  ship  section  is  divided 
into  a  number  of  discrete  straight  line  segments,  as  is  done  for  the  source  distribution  method. 
The  velocity  potential  ^  is  assumed  to  be  constant  along  each  segment.  Equation  (4.1)  is  satisfied 
at  the  midpoint  of  each  segment,  giving  the  following  discretized  solution  of  the  above  equation: 


The  matrix  form  of  Equation  (4.2)  is: 


The  terms  of  Equation  (4.3)  are: 


Gij  —  „  Sij 


[C]  {<!>}  =  {D} 

dG(lxi) 


i-/ 

2^  Js. 


dn^ 


di 


2N 


Di  —  ^  ^  Eij 


i=i 


Eij  — 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


For  symmetrical  ship  sections,  the  symmetry  conditions  of  Equations  (3.15),  (3.16),  (3.17), 
and  (3.18)  apply  and  only  the  N  potentials  for  ?/  >  0  need  to  be  solved.  The  revised  matrix 
coefficients  are: 


^ij 


= 


_  1  Sij  +  — 
2  27r 


_  1  Sii  +  — 

2  27r 


N 


Jsi  Js: 


L 


dG{^,Xi) 


dn? 


+ 


i 


dG{f*,Xi)  J 
dGi{*,Xi) 


5»  dn^ 


EjP  for  A:  =  2,3,4 
j=i 


II 

1 

27r 

^  G{i,xi)d(  -  j 

I - 1 

t* 

II 

1 

27r 

\j^  G{lxi)d(  +  j 

G{f*,xi)df* 

s;  J 

for  A;  =  3 


for  A:  =  2,4  (4.7) 
for  A:  =  3  (4.8) 

(4.9) 

(4.10) 

(4.11) 
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4.1  Improved  Direct  Method 


SclavoTinos  and  Lee  [13]  give  an  improved  direct  method  which  greatly  reduces  irregular 
frequency  problems  of  Equation  (4.1).  Their  method  combines  Equation  (4.1)  with  the  following 
alternative  direct  solution: 


1  d4>{x) 

2  dng 


+ 

27r  drii 


d( 


J_  f  ^  dGjlx)  - 

27r  Js^  dn^  dng 


(4.12) 


The  above  equation  usually  gives  results  that  are  less  accurate  than  Equation  (4.1).  Like  Equa¬ 
tion  (4.1),  Equation  (4.12)  has  irregular  frequencies;  however,  the  two  equations  have  different 
irregular  frequencies. 

The  linear  combination  of  Equation  (4.1)  and  ia  times  Equation  (4.12)  gives  the  following 
solution  for  the  velocity  potentials: 


-  I  + 


2Tr 


d 

dn^ 


^1  +  ia 


Gi(,x)d( 


ia  d(l>{x)  ,  J_  /  d(f>{0  (  d  \ 

2  dng  ^  2-k  4  dn^  \  ^  dng) 


G{ix)d( 


(4.13) 


where  a  is  a  constant  of  ajbitrary  value.  The  solution  of  the  above  equation  satisfies  both  Equa¬ 
tion  (4.1)  and  Equation  (4.12)  if  a  is  non-zero.  An  irregular  frequency  solution  can  satisfy  one 
of  Equation  (4.1)  or  Equation  (4.12),  but  only  a  real  solution  can  satisfy  both  equations  simulta¬ 
neously;  thus,  the  solution  to  Equation  (4.13)  typically  will  not  have  irregular  frequencies  if  a  is 
non-zero.  Sclavounos  and  Lee  recommend  using  a  value  of  0.2  for  a,  which  keeps  Equation  (4.1) 
dominant  but  gives  sufficient  weight  to  Equation  (4.12)  to  mitigate  irregular  frequencies. 

The  discretized  form  of  Equation  (4.13)  is 


-  2  + 


ia  d<j){xi) 
2  dngf 


dm 

dn^ 


^1  -i-  G{(,Xi)d( 


(4.14) 


For  solution  using  Equation  (4.3),  the  revised  matrix  coefficients  (neglecting  symmetry)  are: 


(4.15) 

(4.16) 

(4.17) 
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For  sections  with  lateral  symmetry,  solution  of  velocity  potentials  for  ?/  >  0  can  be  done  using 
the  following  revised  matrix  coefficients: 


rik) 

Di 

Eij 

Eij 


1  A  .  d 

2x  V 

-i/; 

2 


xi  J  Jsj  Js* 


for  A;  =  2,4  (4.18) 


+  ^|l  +  ia 


“5^)  Is, 


dn^ 


+ 


L 


duz 


la 


N 


Vi  +  ^  Vj  Eij  for  A:  =  2, 3, 4 


j=i 


dr  I  for  A:  =  3  (4.19) 

(4.20) 

for  A:  =  2,4  (4.21) 

for  A;  =  3  (4.22) 


For  a  ship  section  with  lateral  symmetry  and  N  segments  on  its  half-section,  the  direct 
method  requires  solution  of  an  iV  x  iV  system  of  equations  regardless  of  whether  Equation  (4.1) 
or  the  modified  form  of  Equation  (4.13)  is  used.  The  main  burden  of  using  the  more  complex 
equation  is  that  additional  derivatives  of  the  Green  function  integrals  are  required.  Fortunately, 
these  additional  derivatives  introduce  relatively  little  new  computational  effort,  as  discussed  in 
the  next  section. 


5  Evaluation  of  Green  Function  and  Its  Derivatives 

The  evaluation  of  the  Green  function  of  Equation  (3.2)  is  a  critical  aspect  of  both  source 
distribution  and  direct  methods.  Among  various  references  discussing  the  two  dimensional  Green 
function  (e.g.  PVank  [7],  Kim  [16],  and  Sarpkaya  and  Isaacson  [17]),  Newman  [19]  gives  the  most 
useful  discussion  for  efficient  implementation  into  a  computer  code. 

The  Green  function  of  Equation  (3.2)  can  be  divided  into  components  which  are  indepen¬ 
dent  and  dependent  of  frequency  as  follows: 

G{x,Xs)  =  G{x,Xs)  +  G{x,Xs)  (5.1) 

G(x,Xs)  =  Inr  —  Inri  (5-2) 

G(x,Xs)  =  —  2  ®  - exp  [— )uZ]  cos{^Y)  djj,  +  i  27r  exp(—Z)  cosY  (5.3) 

Jo  -  1 
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where  G  is  the  frequency  independent  component  and  G  is  the  frequency  dependent  component. 
The  following  two  subsections  discuss  evaluation  of  each  of  the  above  components. 


5.1  Frequency  Independent  Term  of  Green  Function 

The  frequency  independent  component  of  the  Green  function  is  very  easy  to  evaluate; 
however,  the  source  distribution  and  direct  methods  require  integration  of  the  Green  function 
across  a  source  or  field  line.  Experience  has  shown  that  the  frequency  independent  component 
must  be  integrated  accurately  to  obtain  accurate  flow  solutions.  Consequently,  integration  should 
be  performed  analytically  rather  than  using  numerical  schemes  such  as  Gaussian  quadrature. 

To  simplify  the  present  discussion,  only  the  In  r  term  of  the  Green  function  is  considered. 
Treatment  of  the  In  ri  term  is  essentially  the  same. 

The  flrst  integral  to  be  considered  is  that  from  a  line  source  with  end  points  x^i  and  Xs2 
to  a  field  point  x,  as  shown  in  Figure  4.  The  inclination  of  the  source  is: 


6s  =  arctan  ( — — (5.4) 
\ys2  -  VsiJ 

The  angles  of  the  source  end  points  relative  to  the  field  point  are: 


=  arctan  ( 

1  (5.5) 

V 

y-ysij 

Ots2 

=  arctan  ( 

Z-  Zs2^ 

1 

1  (5.6) 

V 

y-ys2j 

|Q!52  -  ttsll 

<  TT 

(5.7) 

According  to  Comeau  and  Ando  [21], 


f£s2  f 

/  In  r  dxs  =  cos  6s  <  {Vsi 
Jx.i  I 


the  integral  across  the  line  source  of  the  In  r  term  is: 

-  ys2)  -  ^{y  -  ys2)ln  [{y  -  ys2f  +  (^  -  2*2)^] 

+  ^(2/  -  2/si)ln  [{y  -  ysif  +  (2  -  2^1  )2]  | 


+  sin 0s  |(2:si  -  Zs2)  -  i(2  -  2s2)ln  [{y  -  ys2f  +  (^  -  ^*2)^] 

+  ^(2  -  2si)ln  [(y  -  ysif  +  (2  -  2si)^]  I 

+  (q:s2  -  ftsi)  [{z  -  Zs2)  cos  0s  -  (j/  -  ?/s2)  sin0s]  (5.8) 
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The  derivatives  of  Equation  (5.8)  with  respect  to  the  field  point  coordinates  are: 

/- 


51nr 

dX  e 


/Xsl 


dy 


cos 


[(y  -  ys2?  +  (2  -  Zs2f]  - 


{y-ys2)^ 


ys2f  +  {z-  Zs2y 


+  |ln[(y-t/si)^  +  (z-2:5i)^]  + 


{<0 


-{z  -  Zs2){y  -  ys2) 


{y  -  ys\f  +  {z  -  ZsxY , 
{z  -  Zsi){y  -  ys\)  \ 


j_  qin  ^  J  V  _ J_  _ 

^  ^  ^  ys2f  +  [z-  z,2f  ^  {y  -  y^  +  {z-  z^ 


+ 


\  dy  dy  ) 


[(2  -  2s2)cos0i  -  {y  -  ys2)sin0s]-  {as2  -  ocsi)  sin 9s  (5.9) 


f^sS 

Jxsl 


^Inr 
dz 


dxs  = 


cos^, 


-{y  -  ys2){z  -  Zst) 


+ 


{y  -  ys\){z  -  Zsi) 

(y-  ysiT  +  {z  -  Zsif  ^ 

(Z  -  Zs2f 


:} 


ys2f  +  (z-  Zs2y 


ys2Y  +  (2:-  2^52)^ 

+  sin(9s|-iln[(y-j/s2)^  +  (2-2i2)^]  - 

~  (^“^*2)sin^s]  +  (as2  -  asi)  cos  Os  (5.10) 


The  derivatives  of  angles  0:51  and  are: 

doisi  _ 
dy 

dosi 


-{Z-Zsx) 


{y  -  ysif  +  {z  -  Zsif 

_  ^  y-y$\ 

dz  {y  -  ysiy  +  {z  -  Zs^y 

dOis2  _  -{z-Zs2) 

dy  ~  {y-  ys2y  +  {z-  Zs2y 

das2  _  y-ys2 _ 

dz  (y  -  ys2y  +  (z-  Zs2f 

The  integral  from  a  point  source  Xs  to  a  field  line  with  end  points  xi,X2  shown  in 
is  also  required.  The  field  line  has  parameters  similar  to  the  line  source  as  follows: 


6  =  arctan 


22-^1 

y2-yi 


(5.11) 

(5.12) 

(5.13) 

(5.14) 
Figure  5 

(5.15) 
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z 


y  - 

Field  point 

•  X 


xli 


Xl2 


Figure  4:  Line  Source  and  Field  Point 


ai 

=  arctan  ( 

^  Zs-  Zl\ 

(5.16) 

\ 

Kys  -  2/1/ 

0^2 

=  arctan  1 

^  Z$-  Z2'\ 

(5.17) 

\ 

^ys  -  2/2/ 

\a2  -  ai\ 

<  TT 

(5.18) 

The  integral  across  the  field  line  of  Inr  is  essentially  the  same  as  Equation  (5.8).  The  integral 
of  the  normal  derivative  of  In  r  along  the  field  line  is  required  for  the  direct  method.  When 
considering  the  orientation  of  the  normal  vector  on  the  field  line,  it  is  assumed  that  field  line 
points  are  numbered  consecutively  outward  from  the  centerline  y  =  Q  and  that  the  normal  vector 
to  the  field  line  points  into  the  ship.  From  Comeau  and  Ando  [21],  the  integral  of  the  normal 
derivative  on  the  field  line  is  simply: 


L 

I. 


ainr  _ 

■fl — 
xi 


d  hi.  r 

X\ 


ai  -  02  for  yi>0  and  j/2  >  0 

—  (oi  -  02)  for  2/1  <  0  and  2/2  <  0 


(5.19) 

(5.20) 


The  derivative  of  Equation  (5.19)  with  respect  to  the  normal  vector  at  the  source  point  is: 


d  51nr  / dai  da2\  ( da\  da2\ 


(5.21) 


The  derivative  of  Equation  (5.20)  can  be  similarly  derived. 

Special  attention  must  be  given  to  evaluation  of  the  In  r  integrals  when  the  point  x  under 
consideration  is  the  midpoint  of  the  line  segment,  leading  to  singularities.  In  this  situation,  the 
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y 


Point  source 


Xi 


^2 


Figure  5:  Field  Line  and  Point  Source 


following  solutions  exist: 


d 


I 


r^s2 

/  In  r 
Jxsl 

^s2  Qlxir 


dn^t 


fXsl 

51nr 

X. 


fXsX 


dng^ 


dxs 


dxg 


dxg 


(5.22) 

(5.23) 

(5.24) 


where  Is  is  source  length.  No  such  singularities  exist  with  the  term  In  ri,  except  at  the  waterline 
2  =  0  where  In  r  and  In  ri  cancel  each  other. 


5.2  Frequency  Dependent  Term  of  Green  Function 

The  frequency  dependent  term  of  the  Green  function  is  more  complicated  than  the  fre¬ 
quency  independent  term.  Fortunately,  its  integration  over  a  line  source  or  field  line  can  be 
accurately  approximated  using  Gaussian  quadrature  with  a  small  number  of  points  on  the  line 
(typically  less  than  eight). 

Following  Newman  [19],  the  frequency  dependent  term  of  the  Green  function  is  evaluated 
as: 

G{x,Xs)  =  -2  Real  {Fpy(y,  Z)}  -|-  i27r  exp(— Z)  cos(y)  (5.25) 

where  Fpv{Y,  Z)  is  a  function  representing  the  principal  value  integral  of  Equation  (5.3).  New¬ 
man  introduces  the  following  complex  variable  Zc'. 

Zc  =  -Z  -I-  iY  (5.26) 


14 


The  function  Fpv{Y,  Z)  can  be  re-written  as: 

FpviY,Z)  =  ex-pZcEi{zc)  +  iriexpzc 
where  Ei{zc)  is  the  exponential  integral  defined  by: 


i(^c)  =  r 

J  Zr. 


dt  ,  I  arg(2c)|  <  tt 


(5.27) 


(5.28) 


Newman  gives  three  efficient  methods  for  evaluating  Equation  (5.27)  based  on  the  ranges  of  Y 
and  Z  shown  in  Figure  6.  When  \zc\  is  small  (T  <  4  and  Z  <  8),  the  following  ascending  series 
gives  rapid  convergence: 


Fpv{Y,Z)  =  expzc  -7  -  In  2c  -  ^ 


(5.29) 


The  above  series  requires  approximately  thirty  terms  to  obtain  convergence  to  six  decimal  places. 

For  large  values  of  Y  or  \zc\  (Y  >  4  or  SY  >  32  -  2Z),  the  exp  Zc  Ei{zc)  term  of  Equa¬ 
tion  (5.27)  is  evaluated  using  a  continued  fraction  expansion: 


expzc  Ei{zc) 


(5.30) 


2c  +  .  .  • 

A  disadvantage  of  the  continued  fraction  expansion  is  that  it  is  not  possible  to  specify  a  tolerance 
and  then  to  compute  successive  terms  until  that  tolerance  is  obtained.  Instead,  the  number  of 
terms  in  the  continued  fraction  must  be  decided  upon  before  computing  terms.  Fortunately, 
excellent  convergence  is  obtained  using  40  terms  for  the  range  of  Y  and  Z  given  in  Figure  6. 

A  double  series  expansion  is  used  for  small  values  of  Y  and  intermediate  values  of  Z  (Y  <4 
and  3y  <  32  —  2Z  and  Z  >  8): 


Fpv(Y,Z)  =  £ 


(iY)-  [^(m-1)! 


—  exp(— .^)  Ei{Z) 


where  Ei{Z)  is  a  second  exponential  integral  from  Abramowitz  and  Stegun  [22]: 


Ei{Z) 


(Z>0) 


(5.31) 


(5.32) 


The  exponential  integral  Ei{Z)  is  easily  evaluated  using  the  following  series  representation  from 
Reference  22: 


Ei{Z)  =  7  -b  InZ  -b  £ 


(^>0) 


(5.33) 
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Figure  6:  Ranges  for  Evaluation  of  Principal  Value  Integral  of  Green  Function 


Evaluation  of  Fpv{Y,  Z)  using  Equation  (5.31)  normally  converges  to  6  decimal  accuracy  within 
n  =  20  terms. 

The  restriction  of  Equation  (5.28)  that  ]  arg(zc)|  <  suggests  that  Y  cannot  be  zero 
(i.e.  the  source  and  field  points  cannot  have  the  same  horizontal  coordinates).  In  practice,  the 
ascending  series  (Equation  (5.29))  and  double  series  (Equation  (5.31))  allow  Fpv{Y,Z)  to  be 
evaluated  for  Y  =  0.  The  only  restriction  on  Fpv(Y,  Z)  is  that  either  Y  or  Z  must  be  non-zero. 
Noting  that  Z  is  zero  only  when  both  the  source  and  field  point  are  on  the  waterline  2  =  0, 
the  problematic  situation  of  both  Y  and  Z  being  zero  only  occurs  when  the  source  point  is 
located  on  the  waterline  and  is  the  same  as  the  field  point.  This  can  occur  when  using  a  deck 
lid  to  remove  irregular  frequencies  in  the  source  distribution  method.  During  computations  the 
problem  can  be  avoided  by  using  Gaussian  quadrature  when  evaluating  the  influence  of  a  line 
source  at  the  waterline  on  the  flow  at  the  line  midpoint. 

The  derivatives  of  the  frequency  dependent  part  of  the  Green  function  are  easily  evaluated 
once  Fpv{Y,  Z)  is  known.  Correcting  a  sign  error  in  [19],  the  derivatives  of  Fpv(Y,  Z)  are: 


dFpv 

~W~ 

dFpv 

dZ 


i  Fpv  + 
-Fpv  - 


i 

Z  -  iY 
1 

Z  -  iY 


Additional  derivatives  of  Fpv{Y,Z)  are: 


d^Fpv 

=  -  Fpv 

1 

1 

dY^ 

Z  -  iY 

(Z  -  iYf 

d^Fpv 

—  Fpv  + 

1 

1 

1 

dZ'^ 

Z  -  iY 

(Z  -  iY)2 

d'^Fpv 

=  -i  Fpv 

i 

i 

dYdZ 

1 

N 

1 

(Z  -  zY)2 

(5.34) 

(5.35) 

(5.36) 

(5.37) 

(5.38) 


Using  the  above  derivatives,  the  derivatives  of  the  frequency  dependent  part  of  the  Green  func¬ 
tion  are: 


dG{x,  Xs) 
dy 


dG{x,  Xs) 
dz 

d‘^G{x,  Xs) 
dydvs 


-2  k  Sign  {y  -  ys) 

X  Real  |eTpv(y,  ^  ^  iY  }  ^  exp(-Z)  sinY 
Real|-f>v(Y,Z)--^^} 


2k 

2P 


+  iir  exp(— Z)  cosY 


(5.39) 

(5.40) 


Real 


|-E>v(Y,  Z)  -  exp(-Z)  cos  Y  (5.41) 
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d^G{x,x,) 

dydzs 


d^Gjx,  Xs) 
dzdys 

d^G{x,  Xs) 

dzdzg 


=  2k'^  Sign{y-ys) 


Real  -iFpv{Y,  Z) 


d^G(x,  Xs) 
dydzs 

d^G(x,Xs) 

dydzs 


i 

Z-iY 


(Z  -  iY)^ 


i-K  exp(— Z)  siny 


(5.42) 


(5.43) 

(5.44) 


Unlike  the  frequency  independent  part  of  the  Green  function,  the  frequency  dependent  part 
and  its  derivatives  can  be  integrated  to  sufficient  accuracy  using  Gaussian  quadrature.  Cook 
[23]  gives  a  good  summary  of  quadrature  techniques.  As  an  example,  the  integral  of  G  from  a 
line  source  can  be  approximated  by: 

Ng 

G{x,  Xs)  dxs  —  l.  Y, 
i=i 

where  No  is  the  number  of  quadrature  points  on  the  line  segment,  Wj  is  the  weighting  factor  of 
point  y ,  and  Xg-j  is  the  location  of  quadrature  point  j  on  the  line  source.  Figure  7  shows  a  line 
source  with  four  quadrature  points.  The  sum  of  the  weights  Wj  is  equal  to  one.  For  solution  of 
hydrodynamic  coefficients  using  the  source  distribution  and  direct  methods,  between  one  and 
four  quadrature  points  per  line  segment  typically  will  give  sufficient  accuracy.  Using  an  even 
number  of  quadrature  points  (JVg  =  2  or  4)  eliminates  the  problem  of  both  Y  and  Z  being  zero 
when  evaluating  the  flow  induced  by  a  deck  lid  source  on  its  own  midpoint. 


z 


y 


Field  point 

•  X 


Figure  7;  Gaussian  Quadrature  Schematic  for  Integration  of  Green  Function 
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6  Added  Mass  and  Velocity  Potentials  at  Low  and  High 
Frequency  Limits 


The  solutions  for  velocity  potentials  and  hydrodynamic  coefficients  presented  thus  far  do 
not  permit  solution  at  the  low  frequency  limit  w  0  or  at  the  high  frequency  limit  u>  oo.  A 
solution  for  a;  =  0  can  provide  a  useful  approximation  for  low  frequency  oscillations.  A  corre¬ 
sponding  solution  oj  =  oc  can  be  useful  for  approximating  high  frequency  oscillations  for  which 
numerical  problems  can  occur  with  the  frequency  dependent  solutions  presented  previously. 

To  solve  for  added  mass  and  velocity  potential  at  the  low  and  high  frequency  limits,  the 
time-independent  velocity  potential  of  Equation  (2.8)  is  replaced  by: 


±  (6.1) 

t  CjJ 

The  body  boundary  condition  of  Equation  (2.11)  simplifies  to: 

=  nW  on  56  (6.2) 

on 

The  solution  of  (f>'  can  be  found  using  either  the  source  distribution  method  or  direct  method 
described  earlier. 

Garrison  and  Berklite  [24],  Garrison  [25],  and  Newman  [26]  discuss  the  solution  of  added 
mass  and  velocity  potentials  at  low  and  high  frequency  linoits.  The  surface  body  condition  of 
Equation  (2.4)  can  be  expressed  in  terms  of  4>'  as: 


-  +  g  =  0  at  z  =  0  (6.3) 

az 

The  low  frequency  limit  of  this  boundary  condition  is: 

=  0  at  z  =  0  for  w  =  0  (6.4) 

oz 

The  high  frequency  limit  of  the  boundary  condition  is: 

4/  =  0  at  z  =  0  for  OJ  =  (XI  (6.5) 


The  Green  function  that  satisfies  the  zero  frequency  boundary  condition  Equation  (6.4)  for  sway 
motion  is: 


Go{S,()  =  lnr-1-lnr-i  (6.6) 

For  heave  and  roll  of  a  surface-piercing  body,  Newman  indicates  that  the  radiation  condition  of 
Equation  (2.5)  is  no  longer  valid  at  zero  frequency.  Unfortunately,  the  Green  function  for  heave 
and  roll  of  a  surface-piercing  body  at  the  zero  frequency  limit  is  unknown. 

The  Green  function  satisfying  the  infinite  frequency  boundary  condition  of  Equation  (6.5) 
is  the  frequency  independent  component  of  Equation  (3.2): 

Goo{S,Xs)  =  Inr-lnri  (6.7) 


19 


Since  the  Green  functions  for  both  zero  and  infinite  frequency  are  real,  the  potential  <!>'  is  real. 
Added  mass  is  easily  computed  using  the  following  revised  form  of  Equation  (2.16): 


^iij 


(6.8) 


Approximations  to  the  complex  velocity  potential  ^  at  low  and  high  frequencies  can  be  made 
using  Equation  (6.1),  with  the  exceptions  of  heave  and  roll  for  surface-piercing  bodies  at  low 
frequencies. 

XJrsell  [27]  discusses  the  heave  motion  at  low  frequencies  of  elliptical  sections  which  are 
half-submerged.  The  added  mass  can  be  approximated  by: 


M33 


pB^ 


(6.9) 


where  B  is  the  waterline  breadth.  Ursell  gives  the  following  equation  for  the  coefficient  K33: 


Kzz 


In 


0.23 


(6.10) 


where  T  is  sectional  draft.  For  other  surface-piercing  sections,  one  can  estimate  Kzz  using 
Equation  (6.9)  if  added  mass  pzz  is  known  at  a  low  frequency.  Equation  (6.9)  can  then  be  used 
to  estimate  heave  added  mass  at  other  low  frequencies 


7  Numerical  Implementation  of  Source  Distribution  and 
Direct  Methods 

Subroutines  have  been  developed  for  computation  of  hydrodynamic  coefficients  using  the 
source  distribution  and  direct  methods.  The  most  difficult  aspect  of  the  numerical  implementa¬ 
tion  was  assembling  the  equations  presented  earlier  in  a  consistent  manner.  Different  references 
use  different  conventions  for  various  terms,  creating  problems  in  consistent  numerical  implemen¬ 
tation.  For  example,  some  references  put  a  factor  of  2/7r  in  front  of  the  Green  function  for  a  two 
dimensional  source.  Sign  conventions  for  axes  and  normal  vectors  introduce  additional  sources 
of  possible  errors. 

After  obtaining  the  relevant  equations  in  a  consistent  manner,  numerical  implementation 
of  the  source  distribution  and  direct  methods  was  relatively  simple.  The  majority  of  developed 
source  code  is  used  by  both  methods;  thus,  the  incremental  effort  required  for  having  both 
methods  available  was  minimal. 

For  integrations  of  Green  functions  across  line  segments,  terms  involving  Inr  and  Inrj 
should  be  integrated  carefully,  particularly  if  t  and  r\  are  small.  The  implemented  code  inte¬ 
grates  terms  with  In  r  and  In  ri  exactly  for  all  cases.  For  the  frequency  dependent  part  of  the 
Green  function,  Gaussian  quadrature  using  four  quadrature  points  on  the  line  segment  gives 
sufficient  accuracy. 


20 


8  Hydrodynamic  Coefficients  for  Sample  Geometries 


This  section  gives  hydrodynamic  coefficients  for  various  sample  geometries.  The  sample 
computations  provide  useful  information  regarding  irregular  frequencies,  differences  between  the 
various  methods,  and  behaviour  at  lower  and  upper  frequency  limits. 


8.1  Sectional  Geometries  for  Computing  Hydrodynamic  Coefficients 

Eight  different  sectional  geometries  are  used  for  testing  the  source  distribution  and  direct 
methods.  The  geometries  shown  in  Figures  8  to  15  represent  a  broad  range  of  sections  for  which 
hydrodynamic  coefficients  are  required.  Table  1  gives  the  offset  coordinates  for  the  sections. 
The  first  three  geometries  are  representative  of  the  bow,  midship,  and  transom  sections  of  the 
frigate  example  in  the  SHIPM06  manual  [6].  The  fourth  section  is  a  barge  designed  to  test 
the  robustness  of  the  algorithms  for  horizontal  and  vertical  line  segments  for  which  the  close- 
fit  method  can  cause  problems.  The  fifth  section  is  a  submarine  cross-section  which  includes 
the  submarine  sail  (validation  of  the  SUBM02  code  in  Reference  28  indicates  that  SUBM02 
overpredicts  roll  excitation  moment  on  a  submarine  with  a  sail).  The  last  three  sections  are 
circles  at  various  depths  of  submergence. 


CG 


Figure  8:  Bow  Section  of  Frigate 


CG 


Figure  9:  Midship  Section  of  Frigate 
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Figure  10:  Transom  Section  of  Frigate 


CG 

© 


Figure  11:  Barge  Section 


Figure  12:  Submarine  Section  with  Sail 
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Table  1:  Sectional  Offsets  for  Hydrodynamic  Coefficients 

Bow 

y  0  1  1.4  1.8  2.2  2.6  2.8 
z  -4  -3.4  -3  -2.2  -1.4  -0.4  0 
Midship 

y  0  2  4  5  5.6  6  6.2  6.4  6.4 
z  -4  -3.8  -3.2  -2.8  -2.2  -1.8  -1  -0.2  0 
Transom 

y  0  1.4  1.8  2.2  2.4  2.8  2.9  3 
z  -1  -0.9  -0.85  -0.8  -0.6  -0.4  -0.2  0 
Barge 

y  02466666 
z  -4  -4  -4  -4  -3  -2  -1  0 
Submarine  with  sail 

y  0  0.06  0.11  0.165  0.1775  0.165  0.11  0.06  0.0416  0.0416  0.0416  0.0416  0.0416  0.03  0.015  0 

z  -0.708  -0.697  -0.669  -0.595  -0.53  -0.465  -0.391  -0.363  -0.357  -0.33  -0.23  -0.13  -0.048  -0.048 
-0.048  -0.0481 

Deep  circle 

y  0.000  1.042  2.052  3.000  3.857  4.596  5.196  5.638  5.909  6.000  5.909  5.638  5.196  4.596  3.857 
3.000  2.052  1.042  0.000 

z  -18.000  -17.909  -17.638  -17.196  -16.596  -15.857  -15.000  -14.052  -13.042  -12.000  -10.958 
-9.948  -9.000  -8.143  -7.404  -6.804  -6.362  -6.091  -6.000 


Shallow  circle 

y  0.000  1.042  2.052  3.000  3.857  4.596  5.196  5.638  5.909  6.000  5.909  5.638  5.196  4.596  3.857 
3.000  2.052  1.042  0.000 

z  -12.000  -11.909  -11.638  -11.196  -10.596  -9.857  -9.000  -8.052  -7.042  -6.000  -4.958  -3.948 
-3.000  -2.143  -1.404  -0.804  -0.362  -0.091  0.000 

Surface-piercing  circle 

y  0.000  1.042  2.052  3.000  3.857  4.596  5.196  5.638  5.909  6.000  5.909  5.638  5.196  4.596  3.857 
3.317 

z  -11.000  -10.909  -10.638  -10.196  -9.596  -8.857  -8.000  -7.052  -6.042  -5.000  -3.958  -2.948  -2.000 
-1.143  -0.404  0.000 
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Figure  15:  Surface-Piercing  Circle  Section 


8.2  Graphs  of  Computed  Hydrodynamic  Coefficients 

Figures  16  to  57  give  added  mass  and  damping  coefficients  for  the  various  sections.  Water 
density  p  is  not  factored  into  the  plotted  values.  Computational  methods  include  the  direct 
method  (a  =  0),  modified  direct  method  (a  =  0.2),  source  distribution  method  (no  lid),  and 
source  distribution  method  with  lid.  The  source  distribution  method  with  lid  is  used  only  for 
surface-piercing  sections,  with  the  lid  represented  by  two  line  segments. 

A  ®  symbol  at  the  lower  and  upper  frequency  limits  gives  zero  and  infinite  frequency  added 
mass  where  applicable.  Each  page  shows  added  mass  and  damping  plotted  to  the  same  scale  for 
a  given  section  and  mode.  The  damping  is  divided  by  frequency  so  that  the  added  mass  and 
damping  figures  for  a  given  section  and  mode  give  the  relative  magnitudes  of  the  added  mass 
and  damping  forces.  Roll  added  mass  and  damping  are  zero  for  the  circular  sections;  thus,  no 
plots  are  shown. 


Frequency  u  (rad/s) 


Figure  20:  Roll  Added  Mass  for  Bow  Section 


- Direct  (a  =  0) 

- Modified  direct  (a  =  0.2) 

. Source  distribution 

.  Source  distribution  with  lid 


Figure  21:  Roll  Damping  for  Bow  Section 
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Source  distribution  with  lid 


Figure  26:  Roll  Added  Mass  for  Midship  Section 


Figure  27:  Roll  Damping  for  Midship  Section 
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Figure  28:  Sway  Added  Mass  for  Transom  Section 
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Figure  33:  Roll  Damping  for  Transoni  Section 
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Figure  38:  Roll  Added  Mass  for  Barge  Section 


400 


300  - 
Roll 

Direct  (a  =  0) 

Modified  direct  (a  =  0.2) 
Source  distribution 
Source  distribution  with  lid 


Frequency  oj  (rad/s) 


Figure  40:  Sway  Added  Mass  for  Submarine  Sail  Section 
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Figure  41:  Sway  Damping  for  Submarine  Sail  Section 
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Figure  42:  Heave  Added  Mass  for  Submarine  Sail  Section 


Figure  43:  Heave  Damping  for  Submarine  Sail  Section 


Figure  45:  Roll  Damping  for  Submarine  Sail  Section 
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Figure  46:  Sway  Added  Mass  for  Deep  Circle  Section 


- Direct  (o:  =  0) 

- Modified  direct  (a  =  0.2) 
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Figure  47:  Sway  Damping  for  Deep  Circle  Section 
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Figure  SOlrSway  Added  Mass  for  Shallow  Circle  Section 
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Figure  51;  Sway  Damping  for  Shallow  Circle  Section 


Figure  53:  Heave  Damping  for  Shallow  Circle  Section 


Figure  56:  Heave  Added  Mass  for  Surface-Piercing  Circle  Section 


Figure  57:  Heave  Damping  for  Surface-Piercing  Circle  Section 


8.3  Discussion  of  Computed  Hydrodynamic  Coefficients 

In  general,  the  computed  hydrodynamic  coefficients  behave  as  expected,  with  good  agree¬ 
ment  between  the  direct  and  source  distribution  methods.  The  modified  direct  method  clearly 
gives  the  most  reliable  results.  The  direct  method  (a  =  0)  gives  generally  good  results,  but 
has  some  problems  with  irregular  frequencies.  The  source  distribution  methods  have  problems 
regardless  of  whether  a  lid  is  used;  the  presence  of  a  lid  shifts  the  irregular  frequencies.  Irregular 
frequencies  are  particularly  severe  for  the  transom  section,  which  has  a  high  breadth  to  draft 
ratio. 

The  shallow  and  surface-piercing  circle  sections  exhibit  negative  sway  added  mass.  The 
consistency  of  the  results  suggests  that  these  negative  added  masses  could  be  a  real  physical 
phenomenon. 

The  zero  and  infinite  frequency  approximations  for  added  mass  show  excellent  agreement 
with  the  low  and  high  frequency  behaviour  of  the  frequency  dependent  added  masses.  As 
oscillation  frequency  approaches  zero,  heave  added  mass  of  surface-piercing  sections  approaches 
infinity.  In  contrast,  Kim  [29]  shows  that  heave  added  mass  for  semi-submerged  ellipsoids 
approach  finite  limits  as  frequency  approaches  zero.  The  conclusion  from  this  discrepancy  is 
that  two  dimensional  calculations  greatly  over-predict  the  heave  added  mass  of  a  ship  section 
at  low  frequency.  The  next  section  proposes  a  solution  to  this  problem. 

Unlike  the  close-fit  method  currently  in  SHIPMO,  the  new  source  distribution  and  direct 
methods  are  able  to  handle  sections  with  horizontal  and  vertical  segments,  as  demonstrated  by 
the  consistent  results  for  the  barge  section. 


9  Heave  and  Roll  Coefficients  at  Low  Encounter  Frequency 


As  discussed  earlier,  a  zero  frequency  Green  function  provides  a  suitable  approximation  for 
sway  coefficients  at  low  encounter  frequencies.  Unfortunately,  this  low  frequency  approximation 
is  not  valid  for  heave  and  roll. 


Figures  20  to  45  indicate  that  roll  added  mass  approaches  a  finite  limit  and  roll  damping 
approaches  zero  as  oscillation  frequency  approaches  zero.  This  consistent  behaviour  suggests  that 
the  following  approximations  can  be  used  for  predicting  roll  velocity  potential  and  associated 
coefficients  at  low  frequency: 


/i44(w)  =  /f44(wi^^) 

A24(w)  =  0 
A44(a;)  =  0 


(9.1) 

(9.2) 

(9.3) 

(9.4) 

(9.5) 


where  is  the  lowest  frequency  for  which  roll  coefficients  have  been  computed.  From  an 
analysis  of  the  computed  hydrodynamic  coefficients  in  Figures  18  to  57,  it  is  suggested  that 
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should  be  0.2  y/gjrmax  or  lower,  where  Vmax  is  the  maximum  distance  from  the  roll  origin  to  a 
sectional  offset  point. 

In  contrast  to  roll  added  mass  and  damping,  heave  added  mass  for  surface-piercing  two 
dimensional  sections  approaches  infinity  and  damping  approaches  a  finite  value  at  low  frequency, 
as  illustrated  in  Figures  18  to  57.  Heave  added  masses  for  semi-submerged  ellipsoids  presented 
by  Kim  [16]  indicate  that  heave  added  mass  approaches  a  finite  zero  frequency  value  and  heave 
damping  approaches  zero  at  low  frequency.  Comparison  of  strip  theory  and  three  dimensional 
computations  for  semi-submerged  ellipsoids  in  Figures  58  to  65  suggests  that  strip  theory  gives 
excessively  high  added  mass  and  damping  values  for  surface-piercing  sections  at  low  oscillation 
frequencies.  The  ellipsoids  have  draft  equal  to  half  the  beam;  thus,  all  sections  are  semi-circular. 
The  ellipsoid  with  LfB  =  8  can  be  considered  a  rough  approximation  of  a  frigate,  while  the 
ellipsoid  i/H  =  4  roughly  approximates  a  coastal  patrol  vessel. 

Strip  theory  should  be  able  to  provide  reliable  heave  coefficient  predictions  for  a  slender 
ellipsoid  (e.g.  L/B  =  8).  Heave  coefficient  predictions  for  low  LjB  bodies,  and  pitch  coefficient 
predictions  for  all  bodies,  are  prone  to  errors  caused  by  three  dimensional  effects.  Examining 
the  heave  added  mass  and  damping  predictions  in  Figures  58  and  59  for  the  slender  ellipsoid,  it 
is  proposed  that  sectional  heave  coefficients  deteriorate  when  the  oscillation  frequency  is  below 

=  Q.2^Jgly.max,  where  is  the  lowest  frequency  for  reliable  heave  coefficients  and  ymax  is 
maximum  half-breadth  of  a  given  section.  The  cut-off  frequency  is  referenced  to  ymax  instead  of 
y  at  the  surface  to  account  for  sections  which  are  narrow  at  the  waterline  but  are  significantly 
wider  immediately  below  the  surface. 

Heave  computations  at  can  likely  provide  improved  modelling  of  low  frequency  be¬ 
haviour  for  an  actual  ship.  The  approximating  equations  are  as  follows: 


=  iImagU(^)(a;f))| 

(9.6) 

M33(w) 

(9.7) 

AssCw) 

=  0 

(9.8) 

Corrected  strip  theory  coefficients  in  Figures  58  to  65  indicate  that  the  low  frequency  heave 
correction  leads  to  improved  strip  theory  values. 


10  Implementation  of  Improved  Methods  into  Strip  Theory 
Program  SHIPMO 

The  modified  direct  method  should  replace  the  close-fit  method  presently  used  in  SHIPMO 
in  order  to  remove  erratic  results  caused  by  irregular  wave  frequencies.  It  would  also  eliminate 
problems  caused  by  vertical  and  horizontal  section  segments. 

Appendix  A  gives  documentation  for  the  subroutine  BOUND2D,  which  is  the  Fortran 
implementation  of  the  direct  method.  Subroutine  BOUND2D  and  various  supporting  routines 
can  replace  subroutine  CLOSEFIT  in  the  SHIPMO  code.  Integration  of  BOUND2D  within  the 
code  should  be  quite  simple.  Velocity  potentials  from  BOUND2D  are  the  complex  conjugate  of 
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Dimensionless  Frequency 


Figure  58:  Heave  Added  Mass  for  Semi- Submerged  Ellipsoid,  LfB  =  8 


Figure  59;  Heave  Damping  for  Semi-Submerged  Ellipsoid,  L/B  —  8 


Figure  60:  Pitch  Added  Mass  for  Semi- Submerged  Ellipsoid,  L/B  =  8 


Figure  61:  Pitch  Damping  for  Semi- Submerged  Ellipsoid,  LIB  =  8 
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Kim  3D 
Strip  theory 

Strip  theory  with  correction 


Figure  62:  Heave  Added  Mass  for  Semi- Submerged  Ellipsoid,  LfB  =  A 


Figure  63:  Heave  Damping  for  Semi-Submerged  Ellipsoid,  L/B  =  4 


5] 


Figure  64;  Pitch  Added  Mass  for  Semi- Submerged  Ellipsoid,  L/B  =  4 


Pitch 

Damping 

A55/(Ai2a;) 


Figure  65:  Pitch  Damping  for  Semi- Submerged  Ellipsoid,  L/B  =  4 
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values  from  CLOSEFIT.  The  present  SHIPMO  code  calculates  the  complex  conjugate  of  velocity 
potentials  after  calling  CLOSEFIT;  this  step  will  not  be  necessary  when  BOUND2D  replaces 
CLOSEFIT  because  the  BOUND2D  formulation  is  consistent  with  the  rest  of  the  SHIPMO 
code. 

At  high  encounter  frequencies,  the  approximation  based  on  the  infinite  frequency  Green  : 

function  should  be  used  to  compute  hydrodynamic  coefficients  for  sway,  heave,  and  roll.  For  sway  ! 

at  low  encounter  frequencies,  SHIPMO  can  use  an  approximation  based  on  the  zero  frequency  i 

Green  function.  For  heave  and  roll  at  low  encounter  frequencies,  SHIPMO  can  reasonably  assume 
that  velocity  potential  is  proportional  to  frequency;  however,  the  reference  low  frequency  limit 
for  heave  must  be  carefully  chosen.  Furthermore,  the  infiuence  of  the  heave  low  frequency  i 

approximation  on  excitation  forces  and  ship  motions  should  be  examined  before  implementation 
into  SHIPMO. 

i 

I 

11  Conclusions  | 

Two  improved  methods  exist  for  predicting  hydrodynamic  coefiicients  of  two  dimensional 
sections.  Of  the  two  methods,  the  modified  direct  method  gives  more  consistent  results  than  the 
source  distribution  method.  The  modified  direct  method  effectively  eliminates  problems  caused 
by  irregular  frequencies. 

Assymptotic  solutions  are  available  for  sway,  heave,  and  roll  at  high  frequencies  and  for 
sway  at  low  frequencies.  RoU  coefficients  at  low  frequencies  are  well  behaved.  For  low  frequency 
heave  of  surface-piercing  sections,  two  dimensional  methods  give  unrealistically  large  values  ! 

of  added  mass  and  damping  for  ship  sections.  A  cut-off  low  frequency  provides  an  improved  I 

estimate  of  added  mass  at  low  frequencies,  and  could  possibly  give  improved  motion  predictions. 

In  the  next  revision  of  the  SHIPMO  strip  theory  code,  the  modified  direct  method  should  j 

replace  the  existing  close-fit  method.  This  revision  to  SHIPMO  will  eliminate  errors  introduced  ' 

by  irregular  frequencies.  Low  and  high  frequency  approximations  should  also  lead  to  improved 
SHIPMO  predictions. 
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Appendix 


A  Subroutine  BOUND2D  for  Computing  Sectional 
Hydrodynamic  Coefficients 

Subroutine  B0UND2D  computes  hydrodynamic  coefficients  for  a  two  dimensional  symme¬ 
trical  ship  section.  The  subroutine  uses  a  modified  direct  method,  which  effectively  eliminates 
problems  caused  by  irregular  frequencies  for  surface-piercing  sections. 

The  Fortran  declaration  for  the  subroutine  is: 

SUBROUTINE  BDUND2D(I0LDGE0M,  NP,  YSECT,  ZSECT,  ZCG,  GRAVITY,  WE, 

*  ALPHAB2D,  AB2,  AB24,  AB3,  AB4,  PHI2,  PHIS,  PHI4) 

REAL  YSECT (NP),  ZSECT(NP) 

COMPLEX  AB2,  AB24,  AB3,  AB4,  PHI2(NP-1),  PHI3(NP-1),  PHI4(NP-1) 

The  units  of  the  input  and  output  variables  are  determined  by  the  value  of  GRAVITY 
(typically  9.806  m/s^  or  32.174  ft/s^).  The  input  variables  are: 

lOLDGEOM  Set  to  1  if  same  geometry  as  previous  call  to  subroutine,  0  otherwise, 

NP  Number  of  section  offset  points, 

YSECT  Vector  of  NP  lateral  offsets  (m  or  ft),  beginning  at  keel  (YSECT(I)  >  0), 

ZSECT  Vector  of  NP  vertical  offsets  (m  or  ft)  relative  to  waterline,  beginning  at  keel 

(ZSECT(I)  <  0), 

ZCG  Height  of  reference  point  above  waterline  (m  or  ft), 

GRAVITY  Gravitational  acceleration  (m/s^  or  ft/s^).  GRAVITY  is  typically  9.806  (met¬ 
ric)  or  32.174  (British). 

WE  Encounter  frequency  (rad/s). 

The  hydrodynamic  coefficients  must  be  multiplied  by  water  density  to  get  added  mass  and 
damping.  The  output  variables  are: 

AB2  Real  component  is  sway  added  mass  (m^  or  ft^)  and  imaginary  component  is 

sway  damping  (m^/s  or  ft^/s), 

AB24  Real  component  is  sway-roll  added  mass  (m®  or  ft®)  and  imaginary  component 

is  sway-roll  damping  (m®/s  or  ft®/s), 

AB3  Real  component  is  heave  added  mass  (m®  or  ft®)  and  imaginary  component  is 

heave  damping  (m®/s  or  ft®/s), 

AB4  Real  component  is  roll  added  mass  (m^  or  ft**)  and  imaginary  component  is  roll 

damping  (m^/s  or  ft^/s), 
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PHI2 

PHI3 

PHI4 


Vector  of  sway  velocity  potentials  on  segments  (m/s  or  ft/s), 
Vector  of  heave  velocity  potentials  on  segments  (m/s  or  ft/s), 
Vector  of  roll  velocity  potentials  on  segments  (m^/s  or  ft^/s). 


55 


References 


1.  R.T.  Schmitke  and  B.W.  Whitten,  “SHIPMO  -  A  FORTRAN  Program  to  Predict  Ship 
Motions  in  Waves,”  DREA  Technical  Memorandum  81/C,  October  1981. 

2.  R.  Graham  and  G.  Miller,  “SHIPM02  -  An  Updated  User’s  Manual  for  the  SHIPMO  Com¬ 
puter  Program  Incorporating  Measured  Sea  Spectra  and  Wave  Loads,”  DREA  Technical 
Memorandum  84/G,  May  1984. 

3.  R.  Graham,  “SHIPM03:  Improved  Viscous  Roll  Damping  Predictions  for  the  SHIPMO 
Computer  Program,”  DREA  Technical  Memorandum  86/212,  May  1986. 

4.  R.  Graham  and  C.  Trudelle,  “SHIPM04  -  An  Updated  User’s  Manual  for  the  SHIPMO 
Computer  Program  Incorporating  an  Extended  Hydrostatics  Capability  and  an  Improved 
Viscous  Roll  Damping  Model,”  DREA  Technical  Communication  87/304,  March  1987. 

5.  K.A.  McTaggart  and  R.  Graham,  “SHIPM05:  An  Updated  User’s  Manual  Incorporating 
New  Wave  Spectra  and  Ship-Referenced  Forces,”  DREA  Technical  Memorandum  92/212, 
April  1992. 

6.  K.A.  McTaggart,  “SHIPM06:  An  Updated  Strip  Theory  Ship  Motion  Program,”  DREA 
Technical  Memorandum  93/213,  October  1993. 

7.  W.  Frank,  “Oscillation  of  Cylinders  on  or  Below  the  Free  Surface  of  Deep  Fluids,”  Report 
2305,  Naval  Ship  Research  and  Development  Center,  October  1967. 

8.  W.  Frank  and  N.  Salvesen,  “The  Frank  Close-Fit  Ship- Motion  Computer  Program,”  Report 
3289,  Naval  Ship  Research  and  Development  Center,  June  1970. 

9.  O.M.  Faltinsen,  “A  Study  of  the  Two-dimensional  Added-mass  and  Damping  Coefficients 
by  the  Frank  Close-fit  Method,”  Report  69-10-S,  Det  Norske  Veritas  Research  Department, 
1969. 

10.  J.W.  Bedel,  “Improved  Method  for  the  Calculation  of  the  Added  Mass  and  Damping  Coeffi¬ 
cients  of  Oscillating  Cylinders,”  Technical  Note  216,  Naval  Ship  Research  and  Development 
Center,  November  1971. 

11.  S.  Ohmatsu,  “On  the  Irregular  Frequencies  in  the  Theory  of  Oscillating  Bodies,”  Papers  of 
Ship  Research  Institute  (48)  (1975). 

12.  S.  Ando,  “AMADIFF  -  A  FORTRAN  Routine  for  Irregular- Frequency-Free  Calculation  of 
Added  Mass  and  Damping  of  Two-Dimensional  Floating  Bodies  by  the  Extended  Boundary 
Element  Method,”  DREA  Technical  Memorandum  95/213,  May  1995. 

13.  P.D.  Sclavounos  and  C.  Lee,  “Topics  on  Boundary  Element  Solutions  of  Wave  Radiation- 
Diffraction  Problems,”  in  Fourth  International  Conference  on  Numerical  Ship  Hydrodynam¬ 
ics  (Washington,  1985). 

14.  G.W.  Johnston  ajid  S.P.  Foster,  “Frequency  Domain  Prediction  of  Submarine  Motion  in 
Waves  (SUBMO)  -  Volume  1,”  DREA  Contractor  Report  CR/91/435,  July  1991. 


57 


15.  G.W.  Johnston  and  A.  Skarsgard,  “Frequency  Domain  Prediction  of  Submarine  Motion 
in  Waves  -  SUBM02  Final  Report,”  DREA  Contractor  Report  CR/93/438,  July  1992. 
Limited  Distribution. 

16.  W.D.  Kim,  “On  the  Harmonic  Oscillations  of  a  Rigid  Body  on  a  Free  Surface,”  Journal  of 
Fluid  Mechanics  21  (1965). 

17.  T.  Sarpkaya  and  M.  Isaacson,  Mechanics  of  Wave  Forces  on  Offshore  Structures,  Van  Nos¬ 
trand  Reinhold,  1981. 

18.  C.C.  Mei,  The  Applied  Dynamics  of  Ocean  Surface  Waves,  John  Wiley  &  Sons,  New  York, 
1983. 

19.  J.N.  Newman,  “The  Evaluation  of  Free-Surface  Green  Functions,”  in  Fourth  International 
Conference  on  Numerical  Ship  Hydrodynamics  (Washington,  1985). 

20.  N.  Salvesen,  E.O.  Tuck,  and  0.  Faltinsen,  “Ship  Motions  and  Sea  Loads,”  Transactions, 
Society  of  Naval  Architects  and  Marine  Engineers  78,  250-287  (1970). 

21.  R.A.  Comeau  and  S.  Ando,  “AMAD  -  A  FORTRAN  Program  to  Calculate  Added  Mass 
and  Damping  of  Two-Dimensional  Floating  Bodies  by  the  Boundary  Element  Method,” 
DREA  Technical  Communication  90/302,  March  1990. 

22.  M.  Abramowitz  and  LA.  Stegun,  Handbook  of  Mathematical  Functions,  Dover  Publications, 
New  York,  1965. 

23.  R.D.  Cook,  Concepts  and  Applications  of  Finite  Element  Analysis,  John  Wiley  &  Sons, 
New  York,  1981. 

24.  C.J.  Garrison  and  R.B.  Berklite,  “Impulsive  Hydrodynamics  of  Submerged  Rigid  Bodies,” 
Journal  of  the  Engineering  Mechanics  Division  99(EM1),  99-120  (1973). 

25.  C.J.  Garrison,  “Hydrodynamic  Loading  of  Large  Offshore  Structures,”  in  Numerical  Meth¬ 
ods  in  Offshore  Engineering,  edited  by  O.C.  Zienkiewicz,  O.C.  Lewis,  and  K.G.  Stagg  (John 
Wiley  &  Sons,  Chichester,  England,  1978),  pp.  87-139. 

26.  J.N.  Newman,  Marine  Hydrodynamics,  MIT  Press,  Cambridge,  Massachusetts,  1977. 

27.  F.  Ursell,  “On  the  Rolling  Motion  of  Cylinders  in  the  Surface  of  a  Fluid,”  Quarterly  Journal 
of  Mechanics  and  Applied  Mathematics  11(3),  335-353  (1949). 

28.  T.  Dickinson  and  K.  McTaggart,  “Validation  of  Wave  Excitation  Force  Prediction  for  the 
Submarine  Motion  Code  SUBM02,”  DREA  Technical  Communication  94/208,  May  1994. 

29.  W.D.  Kim,  “On  a  Free-Floating  Ship  in  Waves,”  Journal  of  Ship  Research  pp.  182-191,200 
(1966). 


58 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


59 


DC003  2/06^87 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 

1 3.  ABSTRACT  (a  brief  and  factual  summary  of  the  document.  It  may  also  appear  elsewhere  In  the  body  of  the  document  itself.  It  is  highly 
desirable  that  the  abstract  of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the 
security  classification  of  the  Information  In  the  paragraph  (unless  the  document  Itself  is  unclassified)  represented  as  (S)»  (C),  (R),  or  (U),  It 
is  not  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is  bilingual). 

This  technical  memorandum  presents  two  improved  boundary  element  methods  for 
predicting  hydrodynamic  coefficients  of  two  dimensional  sections.  The  first  method 
is  a  source  distribution  method,  while  the  second  is  a  direct  method  for  determining 
velocity  potential.  A  major  objective  of  these  improved  methods  is  to  eliminate 
errors  due  to  irregular  frequencies,  which  currently  exist  in  the  close-fit  method  used 
by  the  strip  theory  program  SHIPMO.  Both  new  methods  successfiilly  eliminate 
problems  caused  by  horizontal  and  vertical  line  segments  with  the  close-fit  method. 

To  eliminate  irregular  frequencies  of  surface-piercing  sections,  the  source 
distribution  method  employs  additional  sources  on  a  deck  lid  at  the  waterline.  The 
introduction  of  a  deck  lid  eliminates  some  irregular  frequencies  but  introduces  new 
irregular  frequencies.  Alternatively,  a  modified  direct  method  effectively  eliminates 
all  irregular  frequencies.  Future  strip  theory  implementations  should  incorporate 
this  modified  direct  method.  This  technical  memorandum  also  reviews  low  and 
high  frequency  limits  of  hydrodynamic  coefficients. 


14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize  a  document  and 
could  be  helpful  in  cataloguing  the  document.  They  should  be  selected  so  that  no  security  classification  is  required.  Identifiers,  such  as 
equipment  model  designation,  trade  name,  military  project  code  name,  geographic  location  may  also  be  included.  If  possible  keywords 
should  be  selected  from  a  published  thesaurus,  e.g.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified. 
If  it  not  posstole  to  select  indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title). 

added  mass 
damping 

hydrodynamic  coefficients 
strip  theory 

boundary  integral  method 


60 

UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


DCD032/06/87-M 


